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1.  INTRODUCTION 


1 . 1  Background.  Ply  waviness,  a  manufacturing-induced  anomaly  in  composite  materials,  is  a  subject 
that  has  received  considerable  attention  from  the  composites  community  (Adams  and  Hyer  1992;  Bogetti 
et  al.  1991a,  1991b;  Garala  1989;  Highsmith  et  al.  1992;  Hyer  et  al.  1988;  Joitner  1984;  Telgadas  and 
Hyer  1990).  Ply  waviness  has  been  demonstrated  experimentally,  numerically,  and  analytically  to 
significantly  reduce  the  mechanical  perfonmance  of  composite  laminates  and  structures.  Hyer  et  al.  (1988) 
and  Telgadas  and  Hyer  (1990)  have  combined  analytic  and  numerical  models  to  study  ply  waviness  in 
laminated  composite  cylinders  for  hydrostatic  pressure  load-bearing  applications.  Adams  and  Hyer  (1992) 
and  Garala  (1989)  have  experimentally  demonstrated  the  deleterious  effects  ply  waviness  can  have  on  the 
mechanical  performance  of  composite  laminates  and  structures.  Bogetti  et  al.  (1991a,  1991b)  have 
developed  linear-elastic  analytic  models  to  study  the  influence  of  ply  waviness  on  stiffness  and  stiengtli 
reduction  in  composite  laminates.  The  importance  of  shear  material  nonlinearity  in  wavy  ply  analyses  has 
also  been  recognized  (Adams  and  Hyer  1992).  In  addition,  interlaminar  shear  has  been  identified  as  at 
least  one  dominant  failure  mechanism  associated  with  ply  waviness  (Bogetti  et  al.  1991a,  1991b;  Hyer 
et  al.  1988;  Telgadas  and  Hyer  1990). 

1.2  Problem  Statement  In  this  work,  analytic  models  are  developed  to  predict  the  effective 
three-dimensional  nonlinear  stress/strain  mechanical  behavior  of  two  different  types  of  wavy  ply 
geometries.  These  models  can  be  used  to  quantify  the  synergistic  effects  ply  waviness  can  have  on  the 
mechanical  performance  of  composite  laminates  and  filament-wound  structures. 

The  two  wavy  ply  geometries  under  investigation  include  a  [90A)/90)  configuration  and  a  [90/±p/90] 
configuration.  A  schematic  illustrating  ply  waviness  in  the  [90A)/90]  configuration  is  shown  in  Figure  la, 
where  an  undulating  [0]  ply  is  embedded  in  [90]  plies  having  straight  fibers.  A  typical  [±P]  configuration 
is  presented  in  Figure  lb  to  illustrate  the  nature  of  the  ply  undulation  at  a  tow  or  yam  crossover  region 
in  a  filament-wound  structure.  In  the  full  [90/±P/90j  configuration  (shown  later),  the  [±p]  plies  will  be 
embedded  between  two  [90]  plies.  Both  of  these  wavy  ply  geometries  represent  sublaminate 
configurations  which  arc  commonly  employed  in  filament-wound  structures.  The  models  developed  in 
this  work  portray  realistic  ply-level  microstructural  anomalies  (waviness)  characteristic  of  the  filament 
winding  process. 
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(a)  [90/0/90]  Model 


(b)  [±P]  Model 


Figure  1. 


2,  ANALYSIS 


2.1  Geometry  Definitions  for  the  [90/Q/901  Model.  From  the  wavy  ply  conAgumion  depleted  (n 
Figure  la,  a  unit  cell  is  first  idealized  (see  Figure  2a).  A  half-sine  wave  segment  of  the  undulation  is 
selected  as  a  representative  element.  Ply  waviness  is  assumed  to  exist  in  the  x-dircetlon  only.  The  length 
of  the  unit  cell  is  L,  the  total  thickness  is  h,,  and  the  uniplitude  of  the  undulation  is  A.  11)C  unit  cell 
consists  of  an  undulating  [0]  ply  (ply  2),  with  thickness  h^  lying  between  two  I90J  ply  regions  (plica  I 
and  3).  The  [90]  plies  are  aligned  in  the  y-direction  and  are  assumed  to  have  no  undulation. 


Detailed  mathematical  descriptions  of  the  Individual  ply  surfaces  within  the  unit  cell,  cxpiesaed  as  a 
function  of  x,  are  required.  Referring  to  Figure  2a,  the  following  equations  arc  used  to  define  iltc  ply 
surface  distances  from  the  laminate  midplane  along  the  half-sine  wave  ."hnped  undulation,  where  llte  flrKl 
subscript  denotes  the  ply  number  and  the  second  subscript  <t  or  b)  denotes  the  top  or  bniiom  ply  iiirfacc, 
respectively. 
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(a)  [90/0/90]  Unit  Cell 


(b)  Out-of-Plane  Orientation 
Angle  Definition 


Figure  2.  Unit  cell  descriptions. 


h,b  (X)  -  -  h,/2 

hi,  (x)  -  hje  (x) 

h2,  (X)  ■  hjb  (x) 

hj,  (x)  *  h|/2, 


—  A/2  —  ttf/2  + 


-  A/2  +  hj/2  + 


(1) 


A  iocti  oui-of-planc  undulation  angle,  9,  is  also  needed.  This  local  angle  changes  along  the 
undulflllon  (i.c,,  is  a  function  of  x),  but  is  assumed  constant  within  any  one  "incremental"  segment. 
Kcferrlng  to  Figuic  2b,  9  can  be  defined  as  a  function  of  x  as 


9  (x)  ■  tan”' 


(2) 


wlicrc  dx  Is  the  width  of  an  incremental  scgmerit  along  the  unit  cell  undulation  in  the  x-direction  and  dz 
IK  the  correxponding  incncmcntal  change  in  ply  surface  height 


2,2  Geometry  Definitions  for  the  [90/±B/90]  Model.  The  unit  cell  idealization  of  the  crossover  region 
(iiiciudliig  the  |9(!]  plies)  is  shown  in  Figures  3a  and  3b.  Geometric  parameters  of  the  unit  cell  are 
quunlKlcd  by  llic  winding  angle,  p,  defined  with  respect  to  the  x-direction,  the  amplitude  of  the  ply  (tow 
or  yarn)  undulation  (assumed  to  be  equal  to  hf),  the  total  sublaminate  thickness,  h,,  and  the  unit  cell 
length,  1  .  of  the  sinusoidal  cros.sovcr  region.  An  out-of-planc  rotation  angle,  9,  at  the  crossover  region 
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is  assumed  to  exist  only  in  the  [±p]  plies.  The  [±|J]  plies  are  embedded  between  two  [90]  plies  of  equal 
and  constant  thickness  (190)  which  are  assumed  to  have  no  undulation. 

Descriptions  of  the  individual  ply  surfaces  within  the  unit  cell,  expressed  as  a  function  of  x,  are  also 
required.  Referring  to  Figure  3a,  the  following  distances  from  the  laminate  midplane  to  the  ply  surfaces 
along  the  half-sine  wave-shaped  undulations  are  defined 


hit  W  =  -  hf 

^Ib  ~  (hf  +  Ijq) 


h4b  (X)  -  hf 

h4t  (x)  =  (hf  +  t^o) .  (3) 
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The  local  out-of-plane  undulation  angle,  6,  for  each  of  the  [±P]  plies  is  defined  analogous  to  the  [0]  ply 
in  the  [90/0/90]  model  presented  in  Figure  2b. 

2.3  Laminate  Analogy.  The  analysis  presented  below  is  generic  in  that  it  can  be  api^ied  to  each  of 
the  two  types  of  wavy  ply  geometries  described  above.  The  effective  nonlinear  mechanical  behavior  of 
the  unit  cell  is  assumed  to  be  representative  of  the  response  of  the  entire  wavy  ply  configuration.  The  unit 
cell  is  first  divided  into  discrete  incremental  segments  of  equal  width,  dx,  in  the  direction  of  undulation 
(the  x-direction).  A  systematic  laminated  media  analogy  is  then  applied  to  each  individual  segment  before 
the  effective  response  of  the  entire  unit  cell  is  formulated. 

An  effective  stress/strain  constitutive  relationship  for  each  discrete  segment  within  the  unit  cell  is  first 
computed  according  to  a  three-dimensional  laminated  media  analysis  (Chou  et  al.  1972) 

oj  «  Cjj  Tj  for  (i,  j  «  1,  2,  3,  4,  5,  6).  (4) 

The  effective  Cj  j  stifihess  matrices  in  Equation  4  are  defined  according  to  the  following  relations 


E  v- 

k-1 


C.'^  - 
u 


c  c 

'“13  ''3j  ^ 

c  ^ 

'“33 


E?.i 


C33 


-,  k  v^n 
'33 


En  V 

1-1  — T 

* 


'33 


for  (i,  j  ■  1,  2,  3,  6),  (5) 


Cy  -  Cji  -  0  for  (i  -  i.  2,  3.  6;  j  -  4.  5) 


(6) 


and 


CiJ 


Ek-i 


v''  pk 

Ak  y 


En  Y'"  ^pk  pi  pk  pi 

k«l  2-rl-l  ''"44  '-SS  ’“45  *“54  > 


Ak  Al 


for  (i,  j  -  4,  5) , 


(7) 
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where 


-  C 


k 

44 


(8) 


and  n  is  the  total  number  of  plies  in  the  segment,  is  the  ratio  of  the  thickness  of  the  k***  ply  to  the  total 
unit  cell  thickness  (i.e.,  (h^^^  -  h^)/ht),  and  is  the  lamina  stiffness  matrix  for  the  k***  ply  in  the 
segment,  which  is  defined  in  the  global  laminate  (x-y-z)  coordinate  system  and  therefore  accounts  for  any 
rotations  (in-plane  and/or  out-of-plane)  due  to  ply  waviness.  Tensoral  rotations  are  used  to  compute  the 
Cjj 's  from  the  typical  engineering  constants  of  the  straight  lamina  (i.e.,  E^,  E2.  E3,  Vj2.  V13,  V23,  G^2'  ^13' 
and  G23)  and  conesponding  rotation  angles  between  the  global  lamiitate  (x-y-z)  and  principal  material 
(1-2-3)  coordinate  systems  (Lehknitski  1963). 

The  local  hj  values  of  the  plies  (within  each  segment)  are  defined  by  Equation  1  for  the  [90/0^] 
model  and  Equation  3  for  the  [9(y±P/903  model.  The  C- ’s  for  the  [90]  plies  (in  both  models)  are  constant 
along  the  x-direction.  The  C^’s  for  the  [0]  arid  [i^J  plies,  however,  must  take  into  account  the 
out-of-plane  (9(x))  orientations  which  vary  along  the  unit  cell.  For  example,  segments  in  the  [90/0/90] 
model  are  represented  by  a  three-ply  laminate  having  out-of-platre  rotations  (0(x))  in  the  Qj  matrices  for 
the  [0]  ply,  Cij(x).  Segments  in  the  [90/±p/90]  model  are  represented  by  a  four-ply  laminate  having  both 
in-plane  (±p)  rotations  and  out-of-plane  (6(x))  rotations  in  the  Cjj  mauices  for  the  (tP)  (^es,  Cjj(x),  and 
Cjj(x).  Compliance  matrices,  ,  for  each  dx  segment  widiin  the  unit  cell  are  then  computed  through 
matrix  inversion  of  the  segment  stifhiess  matrices,  C^. 

2.4  Nonlinear  Lamina  Constitutive  Relations.  Material  nonlinearity  is  introduced  in  the  principal 
shear  directions  of  the  lamina  constitutive  relations,  23, 13,  and  12.  All  other  mechanical  properties  are 
assumed  constant  The  Ramberg-Osgood  Equation  is  adopted  here  to  represent  material  nonlinearity 
(Richard  and  Blackloak  1969).  The  three-parameter  equation  relating  stress  explicitly  as  a  continuous 
function  of  strain  is  given  by 

tW- - w-  ® 

(1  +  (G  "  •  yfx^  p)^'*' 
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where  (7*  is  the  initisl  linear  shear  modulus,  is  an  asymptotic  shear  stress  value,  and  p  is  a  shape 
parameter. 


Taking  the  derivative  of  Equation  9  with  respect  to  strain  yields  an  expression  for  the  tangent 
modulus,  G(,  as  a  continuous  function  of  strain 


G*  -y 


dY  (1  +<G® -y/t®)  p)‘**'P 


(10) 


Equation  10  is  a  convenient  expression  for  determining  the  tangent  or  instantaneous  shear  moduli  of  the 
lamina  within  the  wavy  ply  configuration.  By  continuously  monitoring  the  principal  shear  strain  level  of 
each  ply  within  each  segment  during  incremental  loading,  the  instantaneous  moduli  can  be  computed 
straightforwardly. 

2.5  Incremental  Formulation.  The  theoretical  basis  for  the  nonlinear  incremental  formulation  is  the 
"linear”  analysis  outlined  above.  The  approach  is  to  incrementally  apply  loads  (defined  in  terms  of 
average  three-dimensional  stresses)  to  the  unit  cell  until  the  desired  load  level  is  reached.  The  FuU 
Newton-Raphson  Method  is  the  incremental  solution  technique  adopted  in  the  analysis  (Bathe  and  Cimento 
1980).  Solution  convergence  is  based  on  a  force  tolerance  criteria. 

At  the  beginning  of  each  load  step,  the  instantaneous  lamina  properties  for  all  plies  in  each  segment 
of  the  unit  cell  are  computed  according  to  Equation  10.  All  other  lamina  properties  are  assumed  constant. 
The  instantaneous  shear  properties  computed  from  Equation  10  are  based  on  the  current  principal  ply 
strain  levels  at  the  beginning  of  the  load  step.  The  instaittaneous  shear  moduli  will  vary  from  ply  to  ply 
within  a  given  segment  and  will  also  vary  from  segment  to  segment  along  the  unit  cell,  depending  on  the 
local  strain  distributions.  (Dnce  the  tangent  shear  moduli  are  determined,  the  instantaneous  compliance 
matrix,  for  each  of  the  segments  within  the  unit  cell  are  computed,  as  outlined  above.  The  effective 

incremental  strains  for  each  individual  segment,  due  to  an  applied  load  increment,  A'Cj,  are  computed 
according  to 

ATj*'  *  Ac^,  (11) 


7 


where  the  superscript  k  is  used  here  as  an  index  to  denote  the  segment  number.  Note  that  the 
superscript  k  is  not  used  on  the  load  increment  A^,  since  it  is  assumed  constant  for  all  the  segments 
within  the  unit  cell.  For  each  applied  load  increment,  the  effective  strain  increment  for  the  entire  unit  cell 
is  computed  as  the  sum  of  all  the  segment  strains  by 

A  Ei  -  E  Aef .  (12) 

k-l 

where  m  denotes  the  total  number  of  segments  in  the  unit  cell. 

The  incremental  strategy  used  calls  for  the  application  of  equal  attd  successive  load  it>....,.Tients,  A^, 
to  the  unit  cell  and  a  continuous  monitoring  of  the  resulting  strain  increments  within  the  unit  cell.  For 
each  load  step  that  is  sailed  to  the  unit  cell,  a  piecewise  linear  increment  in  unit  cell  strain  is  computed 
according  to  Equations  1 1  and  12.  The  nonlinear  behavior  of  the  wavy  configuration  is  represented  by 
an  effective  (average)  stress-strain  curve  generated  fhun  the  unit  cell  response  to  incremental  loading.  The 
total  average  stress  at  any  point  during  the  iiKtemental  loading  is  defined  as  the  cumulative  sum  of  all 
prior  load  increments  applied  to  the  unit  cell.  Mathematically,  the  total  average  stress  acting  on  the  unit 
cell.  AoJ  (for  example,  after  nls  load  steps),  is  given  by 


nil 


(13) 


where  the  superscript  k  is  used  here  as  an  index  to  denote  the  load  increment.  The  correspotxling  total 
effective  unit  cell  strain,  e  ^ ,  is  defined  as  the  cumulative  sum  of  all  prior  unit  cell  strain  increments  up 
to  the  current  load  level.  For  example,  after  nls  load  increments,  the  total  unit  cell  strain  is  defined  by 


(14) 


The  effective  stress/strain  behavior  of  the  unit  cell  is  therefore  defined  by  the  relationship  generated 

*r  y 

between  Oj  and  e  ^  over  the  prescribed  load  path  history. 
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2.6  Failure  Criteria.  The  maximum  stress  failure  criteria  is  employed  in  the  analysis.  During  the 
incremental  loading,  the  total  principal  stress  levels  within  each  ply  of  each  segment  are  calculated.  When 
the  principal  stress  level  in  a  ply  reaches  the  value  of  its  corresponding  allowable,  as  defined  in  the 
maximum  stress  failure  criteria,  the  modulus  associated  with  that  ply  and  failure  mode  is  set  arbitrarily 
small  for  subsequent  calculations.  Accordingly,  the  ply  retains  its  current  stress  level,  but  is  prevented 
from  carrying  additional  stress  during  subsequent  incremental  loading. 

3.  RESULTS 

3.1  Input  Summary.  Wavy  model  input  parameters  required  in  the  analysis  include  the  material 
properties  of  the  composite,  tlie  unit  cell  geometric  parameters,  and  details  of  the  incremental  load 
strategy.  An  .>S^/PEIIK  graphiieAhermoplastic  composite  material  system  is  used  to  illustrate  the 
findings.  Material  properties  and  strength  allowables  for  AS4/PEKK  are  listed  in  Tables  1  and  2, 
respeaively.  The  Ramberg-Osgood  parameters  used  to  represent  the  nonlinear  shear  response  of  the 
AS4/PEKK  system  in  the  23,  13,  and  12  principal  shear  directions  are  summarized  in  Table  3.  The 
corresponding  shear  stress/strain  curves  are  presented  in  Figure  4.  For  this  aiudysis,  the  asymptotic  stress 
values  are  assumed  to  be  equal  to  their  corresponding  ultimate  strength  allowables. 


Table  1.  Material  Properties  for  AS4  Graphite/PEKK 


El 

(MPa) 

1.17x10*  1 

Ej 

(MPa) 

9,66  X  10*  1 

(MPa) 

9.66  X  1()* 

Vi2 

0.36 

Vl3 

0.36 

^23 

0.36 

The  incremental  strategy  requires  one  to  define  a  load  step  size,  Aoj,  and  the  total  number  of  load 
increments  to  be  applied.  As  with  most  incremental  numerical  techniques,  solution  accuracy  is  improved 
by  decreasing  the  load  step  size.  It  is  noted  that  convergence  studies  were  conducted  for  all  the  examples 
presented  in  this  work  to  ensure  accurate  numerical  results. 
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Table  2.  Strength  Allowables  for  AS4  Graphite/PEKK 


XIT 

(MPa) 

1.97  X  Kl* 

X2T 

(MPa) 

3.59  X  10* 

X3T 

(MPa) 

3.59  X  10* 

XIC 

(MPa) 

1.07  X  !()* 

X2C 

(MPa) 

1.83  X  10^ 

X3C 

(MPa) 

1.83  X  lU^ 

S12 

(MPa) 

1.31  X  10? 

S13 

(MPa) 

7.59  X  10* 

S23 

(MPa) 

3.79  X  10* 

Table  3.  Ramberg-Osgood  Shear  Parameters  for  AS4/PEKK 

^2 

(MPa) 

6.76  X  !()* 

Gi3 

(MPa) 

6.76  X  10* 

G23 

(MPa) 

4.34  X  10* 

^12 

(MPa) 

1.31  X  10* 

(MPa) 

7.59  X  10* 

tO 

^3 

(MPa) 

3.79  X  10 

P12 

2.0 

Pl3 

2.0 

P23 

2.0 

3.2  (90/0/901  Model  Predictions.  The  first  case  study  employs  the  190/0/90]  wavy  model  with  the 
following  geometric  parameters;  L  =  0.127  cm,  hf  =  0.0127  cm,  h,  =  0.0381  cm,  and  A  =  0.0127  cm.  The 
unit  cell  configuration  is  iiKrementally  loaded  in  uniaxial  compression  in  the  x-direction  until  ultimate 
failure.  The  following  sections  discuss  the  findings  of  this  case  study. 
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Figured  Nonlinear  lamina  shear  stiess/strain  response  for  AS4/PEKK. 

3.2.1  Comparison  With  Linear  Analysis  Predictions.  The  effective  x-ditection  unit  cell  stressAtrain 
response  is  first  examined.  A  linear-elastic  version  of  this  example  was  conducted  to  provide  a  baseline 
for  comparison  with  the  nonlinear  predictions.  The  shear  moduli  employed  in  the  linear  case  are  the 
Ramberg-Osgood  initial  shear  moduli  (see  Table  3).  The  linear  and  nonlinear  x-ditection  stress/strain 
response  predictions  of  the  unit  cell  are  presented  in  Figure  5.  Failure  in  the  linear  case  is  reached  when 
any  principal  ply  stress  compotKot,  in  any  segment,  teaches  its  corresponding  strength  allowable  (see 
Table  2).  Failure  in  the  nonlitKar  case  is  reached  when  the  unit  cell  can  no  longer  sustain  additional 
incremental  loading  (load  divergence  is  detected).  As  expected,  the  stress/strain  responses  are  similar  at 
low  levels  of  strain  and  diverge  at  higher  strain  levels,  where  nonlinear  material  softening  occurs.  The 
ultimate  load  levels  ate  markedly  different  (19S  MPa  and  312  MPa  for  the  linear  and  nonlinear  cases, 
respectively).  The  failure  mechanism  in  the  linear  case  is  interlaminar  shear  stress  failure,  O13,  which 
occurs  in  the  segment  of  the  wavy  configuration  with  the  greatest  ply  undulation  (maximum  shear  stress). 
In  the  nonlinear  analysis,  the  interlaminar  shear  ply  stress,  0^3,  does  not  reach  its  ultimate  but  rather 
approaches  it  asymptotically.  Consequently,  the  wavy  configuration  is  permitted  to  bear  load  well  beyond 
the  load  level  predicted  to  cause  failure  in  the  linear  analysis.  As  the  additional  load  is  ai^lied,  the 
nonlinearity  of  the  strain  response  increases,  to  a  point  where  divergence  is  detected  (312  MPa).  The 
following  sections  help  to  identify  the  fundamental  mechanisms  contributing  to  the  overall  nonlinear 
response  of  the  wavy  unit  cell  configuration. 
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Figures.  Nonlinear  stress/strain  response  of  the  [9(MW01  model  under  uniaxial  compression. 

3.2.2  Development  of  Segment  Strain  Distributions.  Under  the  constant  average  stress  assumption, 
the  segment  strains  which  develop  during  incremental  loading  are  expected  to  vaiy  along  the  unit  ceil, 
peakiiig  where  the  local  out-of-plane  [0]  ply  undulation  angle  is  a  maximum,  at  x  «  L  /2.  Hie 
development  of  nonlinear  segment  strain  distributions  within  the  unit  cell  is  demonstrated  in  Figure  6  at 
load  levels  between  69  MPa  and  276  MPa.  Note  that  the  large  strains  (nonlinear)  begin  to  develop  in 
segments  with  greater  ply  undulation  and  that  they  become  increasingly  more  nonlinear  with  increasing 
load.  This  increased  nonlinear  behavior  is  due  to  the  fact  that  as  the  inteilaminar  shear  strains  increase, 
the  corresponding  tangent  moduli  decrease  and  reduce  the  local  compliance  matrices  within  the  segments. 

3.2.3  Influence  of  Undulation  Amplitude  on  the  Nonlinear  Response.  The  influence  of  undulation 
amplitude  on  the  nonlinear  x-direction  stress/strain  response  of  the  [9(V[V90]  model  (with  the  same 
geometric  parametera  and  under  the  same  loading  condition  as  described  above)  was  also  examined.  In 
Figure  7,  the  x-direction  stress/strain  responses  of  wavy  unit  cell  conflgurations  with  undulation  aiuiplitudes 
A  =  O.OOCXX)  cm,  0.0063S  cm,  and  0.01270  cm  are  presented.  Note  that  the  case  for  A  =  0.00000  cm  is 
linear  to  failure  and  is  identical  to  what  would  be  predicted  by  classical  laminate  analysis  (S24.8  MPa). 
This  is  the  expected  result  since  there  are  no  undulation  or  nonlinear  effects  in  this  case.  Comparing  the 
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other  eases,  it  is  shown  that  the  ultimate  unit  cell  strength  decreases  significantly  with  increasing 
undulation  amplitude.  It  is  noted,  however,  that  these  ultimate  strengths  are  significantly  higher  than  those 
predicted  by  the  linear  wavy  analysis  because  of  the  previously  discussed  reasons  (Figure  5). 

3.3  f90/±B^1  Crossover  Model  Predictions.  In  our  second  case  study,  we  present  a  comparison  of 
the  nonlinear  and  linear  analysis  predictions  for  the  {9(Vtp/90]  crossover  model.  In  Figure  8,  the 
nonlinear  and  linear  stress/strain  responses  of  the  (90/±pW]  crossover  model  for  P  =  64®,  49®,  and  20® 
are  presented.  The  linear  and  nonlinear  predictions  of  the  initial  moduli  for  Ute  respective  models  are  in 
good  agreement  This  is  expected  since  at  low  load  levels,  linear  material  behavior  is  dominant.  As  with 
the  [90A1^]  wavy  model,  the  nonlinear  analysis  predicts  significantly  higher  ultimate  strength  values  than 
those  of  the  linear  analysis.  Again,  this  is  due  to  the  fact  that  the  domiiumt  failure  mechanism  of  the 
linear  analysis,  interlaminar  shear  failure,  does  not  constinite  catostrophic  failure  in  the  nonlinear  analysis. 
It  is  noted  that  these  [90/±p/90]  crossover  configurations  exhibit  substanbally  more  nonlinearity  than  the 
[90A)/90]  wavy  configuration.  This  is  due  to  the  fact  that  the  in-plane  shear  loading  of  tlie  [±p]  plies 
contribute  an  additional  source  of  nonlinearity. 


Figure  8.  Nonlinear  stress/strain  response  of  the  [9Q^B/901  cross-over  model  under  uniaxial 
compression. 


14 


4.  CONCLUSIONS 


An  analytic  model,  based  on  a  three-dimensional  laminated  media  analysis,  was  developed  to  predict 
the  effective  nonlinear  laminate  behavior  associated  with  ply  waviness.  An  undulating  [0]  ply  in  a 
[90A)/90]  sublaminate  configuration  and  an  undulating  [ip]  ply  in  a  [90/ip/90]  sublaminate  configuration 
were  two  types  of  ply  waviness  considered.  An  incremental  loading  sjategy  was  employed  wherein 
piecewise  linear  solutions  were  superimposed  to  obtain  the  overall  nonlinear  stress/strain  response  of 
composite  laminates  with  wavy  plies. 

The  influence  of  ply  waviness  with  nonlinear  shear  material  response  on  the  mechanical  performance 
of  composite  laminates  was  studied.  The  nonlinear  analysis  revealed  significant  nonlinearity  in  the 
stress/strain  response  of  both  the  wavy  ply  and  crossover  region  configurations  investigated.  This 
nonlinear  response  was  largely  due  to  severe  nonlinear  straining  in  regions  of  the  configuration  where  ply 
undulation  was  greatest  This  is  brought  about  by  the  decreasing  tangent  shear  moduli  associated  with 
high  shear  straining  within  undulating  plies.  Strain  distributions  with  the  wavy  ply  configuration  showed 
that  segment  strains  in  regions  of  the  configuration  with  large  ply  undulations  strained  significantly  more 
than  segments  with  no  undulation,  which  exhibited  the  usual  linear  response.  At  low  load  levels,  the 
nonlinear  and  linear  overall  stress/strain  responses  for  both  the  wavy  ply  and  crossover  configurations  were 
similar.  The  most  significant  finding  of  the  nonlinear  analysis  was  that  nonlinear  ultimate  strength 
predictions  for  the  wavy  ply  and  crossover  configurations  were  significantly  higher  than  linear  predictions. 
Interiaminar  shear  is  the  dominant  failure  mechanism  in  the  linear  analysis  strength  predictions.  In  the 
nonlinear  analysis,  this  failure  mechanism  does  not  constitute  ultimate  failure.  Loading  is  permitted  to 
continue  as  loads  are  redistributed  to  adjacent  plies.  This  redistribution  of  load  results  in  ultimate  strength 
predictions  for  the  nonlinear  analysis  as  much  as  twice  that  predicted  in  the  linear  analysis.  The  models 
developed  in  this  work  portray  realistic  ply-level  microstructural  waviness  characteristic  of  the  filament 
winding  process  and  therefore  can  be  used  to  quantify  the  synergistic  effects  ply  waviness  can  have  on 
the  mechanical  performance  of  filament-wound  structures. 
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No.  of 

Cosiss  Organization 

11  Director 

Benet  Laboratories 
ATTN:  SMCAR-CCB, 

C.  Kitchens 
J.  Keane 
T.  Allen 
J.  Vasilakis 
G.  Friar 
T.  Simkins 
V.  Montvori 
J.  Wrzochalski 
G.  D’Andrea 
R.  Hasenbein 

SMCAR-CCB-R.  S.  Sopok 
Watcrvliet,  NY  12189 

1  Commander 

ATTN:  SMCWV-QAE-Q,  C.  Howd 
Bldg  44  Watervliet  Arrenal 
WatcrvUet.  NY  12189-4050 

1  Commander 

ATTN:  SMCWV-SPM.  T.  McCloskey 
Bldg  25/3,  Watcrvliet  Arsenal 
WatcrvUet,  NY  12189-4050 

1  Commander 

WatcrvUet  Arsenal 

ATTN:  SMCWV-QA-QS,  K.  Insco 

WatcrvUet,  NY  12189-4050 

1  Commander 

Production  Base  Modernization  Activity 
U.S.  Army  Armament  Research, 
Development,  and  Engineering  Center 
ATTN:  AMSMC-PBM-K 
Picatinny  Arsenal,  NJ  07806-5000 

1  Commander 

U.S.  Army  Belvoir  RD&E  Center 
ATTN:  STRBE-IBC,  C.  Kominos 
Fort  Belvoir,  VA  22060-5606 

1  U.S.  Army  Cold  Regions  Research  and 
Engineering  Laboratory 
ATTN:  P.  Dutta 
72  Lyme  Road 
Hanover,  NH  03755 


No.  of 

Copies  Organization 

1  Director 

U.S.  Army  Research  Laboratory 
ATTN:  AMSRL-WT-L,  D.  Woodbury 
2800  Powder  Mill  Road 
Adelphi,  MD  20783-1145 

3  Commander 

U.S.  Army  Missile  Command 
ATTN:  AMSMI-RD,  W.  McCotkle 
AMSMI T'  3-ST,  P.  Doyle 
AMSMI  R  D-ST-CN,  T.  Vandiver 
Redstone  Arsenal,  AL  35898 

2  U.S.  Army  Research  Office 

Dir.,  Math  St  Computer  Sciences  Div. 
ATTN:  Andrew  Crowson 
J.  Chandra 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

1  U.S.  Army  Research  Office 
Engineering  Sciences  Div 
ATTO:  G.  Anderson 
P.O.  Box  12211 

Research  Triangle  Park,  NC  27709-2211 

2  Project  Manager 
SADARM 

Picatinny  Arsenal,  NJ  07806-5000 

2  Project  Manager 

Taiik  Main  Armament  Systems 
ATTN:  SFAE-AR-TMA, 

COL  Bregard 
C.  Kimker 

Picatinny  Arsenal,  NJ  07806-5000 

3  Project  Manager 

Tank  Main  Armament  Systems 
ATTN:  SFAE-AR-TMA-MD, 

H.  Yuen 
J.  McGieen 
R.  Kowalski 

Picatinny  Arsenal,  NJ  07806-5000 
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No.  of 

No.  of 

Copies 

Organization 

Conies  Organization 

2 

Project  Manager 

2 

Naval  Surface  Warfare  Center 

Tank  Main  Armament  Systems 

Dahlgren  Division 

ATTN:  SFAE-AR-TMA-MS, 

CodeG33 

R.  Joinson 

D.  Guziewicz 

Dahlgren.  VA  224488 

• 

Picatinny  Arsenal,  NJ  07806-5000 

1 

Naval  Research  Laboratory 

Code  6383 

1 

Project  Manager 

ATTN:  LWolock 

Tai^  Main  Armament  Systems 

Washington,  DC  20375-5000 

SFAE-AR-TMA-MP,  V'.  Lang 

Picatinny  Arsenal,  NJ  07806-5000 

1 

Office  of  Naval  Research 

Mech  Div  Code  1 132SM 

2 

PEO- Armaments 

ATTN:  Yapa  Rajapakse 

ATTN:  SFAE-AR-PM, 

Arlington,  VA  22217 

D.  Adams 

T.  McWUliams 

1 

Naval  Ordnance  Station 

Picatinny  Arsenal,  NJ  07806-5000 

Advanced  Systems  Technology  Br. 

ATTN:  D.  Holmes 

1 

PEO-Field  Artillery  Systems 

Code  2011 

ATTN:  SFAE-FAS-PM,  H.  Goldman 
Picatinny  Arsenal,  NJ  07806 

Louisville,  KY  i»0214-5245 

1 

David  Taylor  Research  Center 

4 

Project  Manager,  AFAS 

Ship  Structures  and  Protection 

ATTN:  LTCD.EUis 

I^partment 

G.  delCoco 

ATTN:  J.  Corrado,  Code  1702 

J.  Shields 

B.  Machak 

Bethesda,  MD  20084 

Picatinny  Arsenal,  NJ  07806-5000 

2 

David  Taylor  Research  Center 

ATTN:  R.  Rockwell 

2 

Commander 

W.  Phyillaier 

DARPA 

ATTN:  J.  Kelly 

Bethesda,  MD  20054-5000 

B,  WUcox 

5 

Director 

3701  North  Fairfax  Drive 

Lawrence  Livermore  National 

Arlington,  VA  22203-1714 

Laboratory 

ATTN:  R.  Christensen 

2 

Commander 

S.  deTeresa 

Wright-Patterson  Air  Force  Base 

W.  Feng 

ATTN:  WL/FIV,  A.  Mayer 

F.  Magness 

Dayton.  OH  45433 

M,  Finger 

P.O.  Box  808 

2 

NASA  Langley  Research  Center 

Mail  Stop  266 

Livermore,  CA  94550 

ATTN:  AMSRL-VS, 

1 

Director 

W.  Elber 

Los  Alamos  National  Laboratory 

F.  Barlett,  Jr. 

ATTN:  D.  Rabem 

Hampton,  VA  23681-0001 

MEE-13,  Mail  Stop  J-576 

P.O.  Box  1633 

■ 

Los  Alamos,  NM  87545 
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No.  of 

No.  of 

Conies 

Organization 

CoDies  Organization 

1 

Oak  Ridge  National  Laboratory 

1 

Purdue  University 

ATTN:  R.  M.  Davis 

School  of  Aero  &  Astro 

P.O.  Box  2008 

ATTN:  C.T.  Sun 

- 

Oak  Ridge,  TN  37831-6195 

W.  Lafayette.  IN  47907-1282 

1 

Battelle  PNL 

1 

Stanford  University 

ATTN:  M.  Smith 

Department  of  Aeronautics  and 

P.O.  Box  999 

Aeroballisticss  Durant  Building 

Richland.  WA  99352 

ATTN:  S.  Tsai 

Stanford,  CA  94305 

6 

Director 

Sandia  National  Laboratories 

1 

UCLA 

Applied  Mechanics  Department, 

MANE  Dept,  Engrg  IV 

Division-8241 

ATTN:  H.  Thomas  Hahn 

ATTN:  C.  Robinson 

Los  Angeles,  CA  90024-1597 

G.  Benedetti 

W.  Kawahara 

2 

Univ  of  Dayton  Research  Inst 

K.  Peiano 

ATTN:  Ran  Y.  Kim 

D.  Dawson 

Ajit  K.  Roy 

P.  Nielan 

300  College  Park  Avenue 

P.O.  Box  969 

Livennare,  CA  94550-0096 

Dayton.  OH  45469-0168 

1 

University  of  Dayton 

1 

Diexel  University 

ATTN:  James  M.  Whitney 

ATTN:  Albert  S.D.  Wang 

300  College  Park  Ave 

32nd  and  Chestnut  Streets 

Philadelphia,  PA  19104 

Dayton,  OH  45469-0240 

2 

University  of  Delaware 

2 

North  Carolina  State  University 

Center  for  Composite  Materials 

Civil  Engineering  Department 

ATTN:  J.  Gillespe 

ATTN;  W.  Rasdorf 

M.  Santare 

L.  Spainhour 

201  Spencer  Laboratory 

P.O.  Box  7908 

Raleigh,  NC  27696-7908 

Newa^,  DE  19716 

1 

University  oi  Illinois  at 

1 

Pennsylvania  State  University 

Urbana-Champaign 

ATTN:  David  W.  Jensen 

National  Center  for  Composite 

223-N  Hammond 

Materials  Research 

University  Park,  PA  16802 

216  Talbot  Laboratory 

ATTN:  J.  Economy 

1 

Pennsylvania  State  University 

104  South  Wright  Street 

ATTN:  Richard  McNitt 

227  Hammond  Bldg 

Urbana,  IL  61801 

• 

University  Park,  PA  16802 

1 

University  of  Kentucky 

ATTN:  Lynn  Penn 

1 

Pennsylvania  Sute  University 

763  Anderson  Hall 

ATTN:  Renata  S.  Engel 

245  Hammond  Building 

University  Park,  PA  16801 
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Lexington.  KY  40506-0046 

No.  of 

No.  of 

Copies  QnfflniatiPIl 

Cooies  Oreanization 

1 

The  University  of  Texas  at  Austin 

1 

Applied  Composites 

Center  for  Electromechanics 

ATTN:  W.  Griscla 

ATTN:  J.  Price 

333  North  Sixth  Street 

10100  Burnet  Road 

St.  Charles.  IL  6Cd74 

Austin,  TX  787584497 

1 

Brunswick  Defense 

1 

University  of  Utah 

ATTN:  T.  Hanis 

Department  of  Mechanical  and 

Suite  410 

Industrial  Engineering 

1745  Jefferson  Davis  Hwy. 

ATTN:  S.  Swanson 

Salt  Lake  City,  UT  84112 

Arlington,  VA  22202 

1 

Chamberlain  Manufacturing 

2 

Virginia  Polytechnical  Institute 

Corporation 

&  State  University 

Research  and  Development  Division 

Dept  of  ESM 

ATTN:  T.  Lynch 

ATTN:  Michael  W.  Hyer 

P.O.  Box  2335 

Kenneth  L.  Reifsnider 

550  Esther  Street 

Blacksburg,  VA  24061-0219 

Waterloo,  lA  50704 

1 

ARMTEC  Defence  Products 

1 

Chamberlain  Manufacturing 

ATTN:  Steve  Dyer 

Corpora'don 

85-901  Avenue  53 

Research  and  Development  Division 

P.O.  Box  848 

ATTN:  M.  Townsend 

Coachella.  CA  92236 

P.O.  Box  2545 

550  Estllier  Street 

2 

Advanced  Composite  Materials 

Waterloo.  lA  50704 

Corporation 

ATTN:  P.  Hood 

1 

Custom  Analytical  Engineering 

J.  Rhodes 

Systems,  Inc. 

1525  S.  Buncombe  Road 

ATTN:  A.  Alexander 

Greer,  SC  29651-9208 

Star  Route,  Box  4A 

Binistcne,  MD  21530 

‘•j 

Alliant  Techsystems,  Inc. 

ATTN:  J.  B^e 

1 

General  Dynamics 

C.  Candland 

Land  Systems  Division 

L.  Osgood 

ATTN;  D.  Bartle 

R.  Buretia 

P.O.  Box  1901 

R.  Becker 

M.  Swensoi) 

Warren,  MI  48090 

600  Second  St.  NE 

3 

Hercules,  Incorporated 

Hopkins,  MN  55343 

ATTN:  R.  Boc 

F.  Policelli 

1 

Amoco  Performance  Products,  Inc. 

J.  Poesch 

ATTN:  M.  Michno,  Jr. 

P.O.  Box  98 

4500  McGinnis  Ferry  Road 

Alpharetta.  GA  30202-3944 

24 

M'lgna,  UT  84044 

No.  of 

No.  of 

Copies 

Organization 

Copies  Organizadon  I 

3 

Hercules,  Incorporated 

2 

Martin-Marietta  Coip. 

ATTN:  G.  Kuebeler 

ATTN:  P.  Dewar 

J.  Vermeychuk 

L.  Sponar 

B.  Manderville,  Jr. 

230  East  Goddard  Blvd. 

Hercules  Plaza 

Wilmington,  DE  19894 

King  of  Prussia,  PA  19406 

2 

Olin  Corporation 

1 

Hexcel 

Flincitbaugh  Division 

ATTN:  M.  Sheiendich 

ATTN:  E.  Steiner 

IISSS  Dublin  Blvd. 

B.  Stewart 

P.O.  Box  2312 

P.O.  Box  127 

DubUn,  CA  94S68-070S 

Red  Lion.  PA  17356 

1 

lAP  Research,  Inc. 

1 

Olin  Corporation 

ATTN:  A.  Challita 

ATTN:  L.  Whitmore 

2763  Culver  Avenue 

10101  9at  Sl,  North 

Dayton,  Ohio  45429 

St.  Petersburg,  FL  33702 

2 

Institute  for  Advanced  Technology 

1 

Rennsaeler  Polytechnic  Institute 

ATTN:  T.  Kiehne 

ATTN:  R.  B.  Pipes 

H.Fair 

P.  Sullivan 

Troy,  NY  12180 

4030-2  W.  Braker  Lane 

1 

SPARTA,  Inc. 

Austin,  TX  78759 

ATTN:  J.  Glatz 

9455  Towne  Centre  Drive 

1 

Integrated  Composite  Technologies 

San  Diego,  CA  92121-1964 

ATTN:  H.  Perkinson,  Jr. 

P.O.  Box  397 

2 

United  Defense  LP 

Yoik  New  Salem,  PA  17371-0397 

ATTN:  P.  Para 

G.  Thomas 

1 

Interferometries,  Inc. 

1107  Coleman  Avenue,  Box  367 

ATTN:  R.  Latriva,  Vice  President 

8150  Leesburg  Pike 

Vienna,  VA  22100 

San  Jose,  CA  95103 

2 

Kaman  Sciences  Corporation 

ATTN:  D.  Elder 

T.  Hayden 

P.O.  Box  7463 

Colorado  Springs,  CO  80933 

3 

LORAL/Vought  Systems 

ATTN;  G.  Jackson 

t 

K.  Cook 

L. L.  Hadden 

1701  West  Marshall  Drive 

• 

Grand  Prairie,  TX  75051 
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No.  of 

No.  of 

Conies  Oreanization 

Copies  Organization 

Aberdeen  Proving  Ground 

AMSRL-WT-PD(ALQ 

A.  Abrahamian 

64  Dir.  USARL 

K.  Barnes 

ATTN:  AMSRL-CI, 

M.  Berman 

C.  Mermegan  (394) 

H.  Davison 

AMSRL-CI-C, 

A.  Frydman 

W.  Stuiek  (1121) 

T.  U 

AMSRL-CI-CB, 

W.  McIntosh 

R.  Kaste  (394) 

E.  Szymanski 

AMSRL-CI-S, 

AMSRL-WT-T, 

A.  Mart:  (309) 

W.  Morrison  (309) 

AMSRL-SL-B, 

AMSRL-WT-TA, 

P.  Dieu  (328) 

W.  GUlich  (390) 

AMSRL-SL-BA. 

W.  Bmchey  (390) 

J.  Walbcit  (1065) 

AMSRL-WT-TC, 

AMSRL-SL-BL, 

K.  Kimsey  (309) 

D.  Bely  (328) 

R.  Coates  (309) 

AMSRL-SL-I, 

W.  de  Rosset  (309) 

D.  HaskiU  (1065) 

AMSPL-WT-TD. 

AMSRL-WT-P, 

D.  Dietrich  (309) 

A.  Horst  (390A) 

G.  Randers-Pehrm  (309) 

AMSRL-WT-PA, 

J.  HufTington  (309) 

T.  Minor  (390) 

A.  Das  Gupta  (309) 

C.  Leveritt  (390) 

J.  Santiago  (309) 

D.  Kooker  (390A) 

AMSRL.WT-W, 

AMSRL.WT.PB, 

C.  Murphy  (120) 

E.  Schmidt  (120) 

AMSRL-WT-WA, 

P.  PlosUns  (120) 

H.  Rogers  (394) 

AMSRL-WT-PC, 

B.  Moore  (394) 

R.  Fifer  (390A) 

A.  Baran  (394) 

AMSRL-WT-PD. 

AMSRL-WT-WB, 

B.  Bums  (390) 

F.  Brandon  (120) 

W.  Drysdale  (390) 

W.  D’Amico  (120) 

K.  Bannister  (390) 

AMSRL-WT-WC, 

T.  Bogctti  (3M) 

J.  Rocchio  (120) 

J.  Bender  (390) 

AMSRL-WT-WD, 

R.  Murray  (390) 

A.  Zielinski  (120) 

R.  Kirkendall  (390) 

J.  Powell  (120) 

T.  Erlinc  (390) 

AMSRL-WT-WE, 

D.  Hq;>kins  (390) 

J.  Temperley  (120) 

S.  Wilkerson  (390) 

C.  McCall  (390) 

D.  Henry  (390) 

J.  Ibomas  (394) 

R.  Kaste  (390) 

L.  Burton  (390) 

J.  'Rteng  (390) 

•* 

R.  Lieb  (390) 

G.  Gazonas  (390) 

M.  Leadore  (390) 

C.  Hoppel  (390) 
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USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratofy  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your 
comments/answers  to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-585 _ DateofReport  .  October  1994 _ 

2.  Date  Report  Received _ ; _ 

3.  Does  this  report  satisfy  a  need?  (Comment  on  purpose,  related  project,  or  other  area  of  interest  for 

which  the  report  will  be  used.) _ 


4.  Specifically,  how  is  the  report  being  used?  (Infonnation  source,  design  data,  procedure,  source  of 
ideas,  etc.)  _ _ 


S.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved, 
operating  costs  avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ ___ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate 
changes  to  organization,  techxucal  content,  format,  etc.) _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address 
above  and  the  Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


department  of  the  Army 


OFFICIAL  BUSINESS 


BUSINESS  REPLY  MAIL 

FIRST  CLASS  PERMIT  NO  0001,  APG,  MD 

PosttgB  will  IM  paid  by  adcbasasB 

Director 

U.S.  Army  Research  Laboratory 

ATTN:  AMSRL-OP-AP-L 

Aberdeen  Proving  Ground,  MD  21005-5066 


